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Abstract 

A detailed study of the dispersion law of surface excitations in liquid ^He 

at zero temperature is presented, with special emphasis to the short wave 

length region. The hybridization mechanism between surface and bulk modes 

is discussed on a general basis, investigating the scattering of slow rotons from 

the surface. An accurate density functional, accounting for backflow effects, 

is then used to determine the dispersion of both bulk and surface excitations. 

The numerical results are close to the experimental data obtained on thick 

films and explicitly reveal the occurrence of important hybridization effects 

between ripplons and rotons. 
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I. INTRODUCTION 



Superfluid helium is a unique system where undamped excitations can propagate also 
in the short wave length regime. This behaviour characterizes the dynamics of the bulk 
as well as of the surface, and is at the origin of several interesting phenomena. Hybridiza- 
tion between ripplons and rotons and quantum evaporation are important examples of such 
phenomena where the surface of helium is explicitly involved. These peculiar properties of 
liquid helium at low temperature have been the object of several experimental investiga- 
tions in the last years. First results for the dispersion of short wave length ripplons have 
recently become available via neutron scattering experiments on helium films |l[ . Quantum 
evaporation, i.e., the emission of an atom produced by a roton propagating balistically in 
the medium and impinging on the surface, has been also systematically investigated (see for 
example Ref. @). 

The above phenomena have received attention also from a teoretical point of view. How- 
ever, due to intrinsic difficulties associated with a reliable description of the statics and 
dynamics of strongly interacting inhomogeneous systems, the predictions have been so far 
only semi-quantitative. The main reason lies in the fact that both microscopic approaches, 
based on correlated basis functions, and phenomelogical methods, mainly based on density 
functional theory, have been used in the framework of the so called Feynman approximation, 
or equivalent schemes. This approximation significantly overestimates the roton energy and 
this clearly rules out the possibility of a quantitative description of quantum evaporation as 
well as of the dynamics of energetic ripplons, since in both cases the relevant energies are of 
the order of the roton gap. 

In order to obtain a better description it is necessary to include backflow effects in the 
theory. Here, following the spirit of density functional theory, we account for such effects 
through a new phenomenological current dependent interaction || . In the present work we 
use the new functional to provide a quantitative description of the ripplon dispersion, with 
special emphasis to the region near the roton threshold. The study of quantum evaporation 
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will be the object of a future work. 

The work is organized in two parts: in Sec. [Tl] we discuss in details the mechanism of 
the hybridization between ripplons and rotons. In particular we show that this mechanism 
is governed by the properties of reflection of rotons from the free surface. In Sec. |T| we 
present numerical results for the ripplon dispersion obtained by solving the equations of 
time dependent density functional theory with the inclusion of backflow terms. 

II. HYBRIDIZATION BETWEEN RIPPLONS AND ROTONS 

The hybridization between ripplons and rotons consists of a coupling between surface 
and bulk excitations that takes place when the energy of ripplons is close to the one of 
rotons. 

The physical origin of such a phenomenon is clear. Let us consider a surface excitation 
with energy larger than the roton gap A. This surface mode can naturally decay into a 
roton if its momentum parallel to the surface is lower than the roton momentum. Indeed 
this process is possible because only the momentum parallel to the surface is conserved, 
while the orthogonal momentum can be transferred to the surface. 

This decay mechanism provides a natural threshold for the energy of ripplons, given by 
the roton gap A. The occurrence of this threshold has been confirmed by recent experimen- 
tal data with neutron scattering [l] and is clearly expected to influence the dispersion of 
energetic ripplons. Such a phenomenon was first explicitly explored in the case of helium in 
Ref. 0] using the formalism of Green's functions. A similar phenomenon was discussed in 
Ref. || for the case of dieletrics where it results in a damping of surface waves. 

The problem of the interaction between ripplons and bulk excitations in superfluid 4 He 
was investigated several years ago by Edwards et al. f§,[/]]. They wrote the solution for the 
velocity potential of the ripplon in the form 

fk — j{ e i(k x x+k z z-ujt) / j\ 
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where x and z are parallel and orthogonal to the surface, respectively; the wave vector k z is 
a function of u; and k x , defined by the equation 

u = <oi+fi) 1/2 ) (2) 

where e(k) is the dispersion law for bulk excitations. The authors of Ref. |||7j took only 
imaginary solutions of Eq. (|2]) for k z . This corresponds to consider purely damped waves in 



the ^-direction with damping length |fc z | _1 . The solution is then given by k z = iyk^.— k 2 {uj) 
where k(u) is the momentum of bulk excitations with energy uj. Of course this solution 
exists only for k x > k(u). However this is not the only possibility, since Eq. (|2|) can have 
also complex solutions for k z . In particular, when u is below and close to the roton gap A, 
one can use the well known approximation M 

«-A + tt^£ (3, 
to the roton dispersion, yielding the following solution for k z : 



■ Po 



k z ~pi + i^2ti(A-w) (4) 
with pi = \Jp% — fef. This solution corresponds to a penetration depth 

t= . Pl (5) 
Po y/2fi(A - u) 

that becomes larger and larger as uj approaches A at fixed k x . 

To investigate the behaviour of the ripplon dispersion near the threshold A one must learn 
more about the dynamics of slow rotons. For this we take the familiar roton Hamiltonian 

Wo 

yielding dispersion (§) for k close to pq. At given k x Eq. (^) can be written in the form 

H = A + t3f4£ . (7) 

with pi defined above. For the description of the threshold behaviour of ripplons only rotons 
with energy close to A are important. In this case the roton wave function can be written 

as 



*r(z) = h(z) exp[-zpiz] + / 2 (z) exp[zpiz] , (8) 

where and ^(z) are slowly varying function of z. By neglecting higher derivatives of 

fi in the Schrodinger equation one finds that the functions fa are solutions of 

H f fi = ufi , i = l,2 (9) 

with 

H ' = A ~^& < 10 > 

and m* = fipl/pf @. 

In the same scheme the surface mode, with energy close to A, appears as a surface bound 
state of rotons, whose wave function can be written as 

^ r (z) = 2Aexp[-Kz] cosQoiz) (11) 

where 

k = [2m*(A , (12) 

and we have made an explicit choice for the phase. Equation (|iT|) corresponds to taking 

h(z) = f 2 (z) = Aexp[-Kz] . (13) 

Clearly both the scattering and the bound solutions (§) and flTT|) hold only far from the 
surface. The quantity k is a function of k x , which can be found, in principle, by solving the 
proper Schrodinger equation in the surface region. Close to the surface, where \z\ is of the 
order of the surface thickness a, the Hamiltonian (|T0|) has to be changed in order to account 
for the interaction of rotons with the surface. However for slow rotons the solutions f\ of the 
Schrodinger equation vary slowly and the role of the surface can be safely replaced by the 
imposition of suitable boundary conditions at z — 0. The general form of these boundary 
conditions is 
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dz 

dfo 

dz 



2=0 



-an/i(0) - ai 2 / 2 (0) (14a) 



= -«2i/i(0)-a 22 / 2 (0). (14b) 

2=0 



The matrix otij must be hermitian in order to ensure real values for k in the case of bound 
states. Furthermore, due to time reversal invariance, one must have an = a 22 and the 
phases of the functions fi(z) and /^(z) can be always chosen such that a± 2 = a^i i s re& l- 
The quantities an and ai 2 , fixing the boundary conditions at z — 0, are entirely determined 
by the properties of the surface. 

Let us first apply Eq. ([TJ|) to the surface mode for uj < A. In this case the functions fi 
are given by expression flT3| ) so that Eq. ( jTl ) yields 

k = yj2m*(A — uj) = a(oo, k x ) , (15) 

where a(u>, k x ) = an + au. This equation has a solution describing a surface mode if a > 0. 
This scheme provides a natural description of the dispersion near the threshold. Let us 
suppose, for example, that a(A, k x ) = at some point k x = k t < p . Near this point the 
function a can be expanded as a (A, k x ) ~ (3{k x — k t ). Thus the dispersion curve near the 
threshold has a parabolic form: 

A-u;=^(k x -k t y , (16) 

approaching the point kt with a horizontal tangent. 



When uj > A the boundary conditions (0) describe the reflection of rotons at the surface. 
In this case the functions f\ take the form 

fi = exp[— ik z z] + Ri(u) exp[ik z z] (17) 

f 2 = R 2 (uj) exp[ik z z] . (18) 

By definition the coefficients R\ and R 2 are the amplitude for change-mode reflection and 
normal reflection, respectively. After imposing the boundary condition (|T4l) in the threshold 
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region to ~ A and near k t where a <C an, a±2, we find the following results for the amplitudes 
Ri and R 2 : 

Ri = —t^— (19) 

ik z + a 

and 

2 ik z + a ^ ^ 

Note that if k z <C a one has 

. (21) 

a 

This equation shows that when approaches A, and consequently k z — > 0, the change-mode 
reflection amplitude approaches the value —1, while the normal reflection amplitude tends 
to zero |1(J . In this limit the solution (§) of the Schrodinger equation, far from the surface, 



takes the form 

%r(z) = exp[-z(pi + k z )z] - exp[-z(pi - fc^z] (22) 

and corresponds to the reflection of a roton R + , with momentum (k x , 0, p\ + into a roton 
R~, with momentum (k x ,0,pi — k z ) (see Fig. [I]). Notice that the roton R~ has negative 
velocity and, consequently, the roton scatters backwards not only in the z direction, but also 
in the x one. The occurrence of the change-mode scattering, peculiar of rotons, will be 
explicitly confirmed in Sec. |T| through the numerical solution of the equations of the time 
dependent density functional theory. 

The theory developed above is very general. It is reasonable to assume however the 
interaction between ripplons and rotons in Tie , producing the hybridization, to be weak. In 
this case one can proceed further, looking for a theory which describes the crossover region 
between an unperturbed ripplon dispersion, entirely characterized by surface properties in 
the absence of interaction with bulk modes, and the dispersion near the threshold, dominated 
by the hybridization mechanism. Under the assumption of weak interaction one expects (see 
the Appendix) that the function a has the form 



a(co,k x ) = b(k x ) + —^ (23) 

e (k x ) - to 

when uj is close to A. In Eq. (|23"D eo(k x ) is the dispersion law of unperturbed ripplons and 
g is a coupling constant. By using Eq. (^) in Eq. (p~5f ) the dispersion of the surface modes 
becomes 

2 

j2m*(A-u) = b{k x ) + — i- . (24) 

Note that the coefficient g 2 must be positive to ensure a stable solution of Eq. (p4[) when 
the unperturbed ripplons have a damping. Different scenarios can now take place depending 
on the value of the coefficient b: 

• 1) b < 0. In this case the spectrum reaches the u = A threshold at a point k x = k t , as 



assumed in the derivation of result (16). The value of k t is fixed by the equation 

Try a = -Kkt) , (25) 
e (fct) - A 

and the coefficient (3 in Eq. fllBD can be directly calculated from Eq. (p4]) . Clearly the 
point fct lies at the right of the unperturbed ripplon branch. 



2) 6 > 0. Equation fl2~4|) has now a solution also for g = 0. A surface branch exists 
independently of the ripplon dispersion. This branch corresponds to a bound state of 
rotons near the surface. This scenario describes a hybridization between two surface 
modes, ripplons and surface rotons respectively, near the threshold for bulk modes. 

3) If the coefficient b satisfies the condition |6| (m*) 1 / 3 ^ 2 / 3 , then its role can be ne- 
glegted in Eq. (p4|). One obtains that the ripplon dispersion approaches asymptotically 
the roton threshold with the law 

A - w oc - - 9 — . (26) 

(e (k x ) - A) 2 v ; 

This hybridization mechanism was previously discussed in Ref. |3J]. Also note that a 
transition from b small and positive to b small and negative corresponds to a transition 
from a bound surface state to a "virtual level" (see §133). 
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In the next Section we will carry out a full calculation of the ripplon dispersion using a 
density functional approach. The numerical results will be shown to follow the behaviour 
predicted by the first scenario (Eqs. (|T6| ) and (p5|)). 



III. PREDICTIONS OF DENSITY FUNCTIONAL THEORY 



In this section we present the results of a calculation of the dynamics of inhomogeneous 
liquid helium obtained using time dependent density functional theory. In the density func- 
tional approach the internal energy is assumed to be a functional of the one-body densities 
in the form 



E 



dr 8[^*^} 



(27) 



where the wave function for a Bose system at zero temperature, is written as 



tf(r,t) = $(r,t)exp(iS(r,t)) 



(28) 



The real function $ is related to the particle density by p = $ 2 , while the phase S is related 
to the velocity v of the fluid by 



j = pv = — (tfVtf* - r*W) = — VS , 
2m m 



(29) 



where j is the current density, and m is the mass of the ^le atoms. 

The formalism of time dependent density functional has been already applied to describe 
the surface excitations of superfluid helium |l2|-p^|]. In Ref. |13| the effective Hamiltonian 
has been derived microscopically starting from the interatomic potential in the framework 



of the hypernetted chain scheme. Vice-versa the functionals employed in Refs. fl4Jl5|l are of 
phenomenological type. 

The equations of motion of time dependent density functional theory can be derived 
starting from the least action principle 



/•*2 

5 dt I dr 

Jti 







(30) 



making variations with respect to $ and S [[OJ, or, equivalently, with respect to ty*. In the 



latter case one finds a Schrodinger-like equation of the form 

H^ = l ^f , (31) 
where H = 5S/5^J* is an effective Hamiltonian. If one looks for linearized solutions 

*(r,t) = * (r,t)+£*(r,*) (32) 
the Hamiltonian H takes the form 

H = H + 5H . (33) 



The static Hamiltonian H is fixed by the equilibrium state ^ (r,t) = e~ llM4t y p(z) (/i 4 is 
the chemical potential) through the equation -f/o^o = Pa^o- The term 5H is linear in 5^ 
and accounts for changes in the Hamiltonian induced by the collective motion of the system. 
Since H depends explicitly on the wave function \I>, the Schrodinger equation ( |3~1~D has to be 
solved using a self-consistent procedure, even in the linear limit considered in the present 
work. The formalism then concides with the one of the Random Phase Approximation 
(RPA) for Bose systems. Of course this theory, which is basically a mean field theory, 
accounts for collective and single particle excitations, but cannot account for multiphonons 
excitations. 

In the bulk \l/ is constant and 5^ can be expanded in plane waves, corresponding to the 
propagation of the phonon-roton mode. In the presence of a free surface, orthogonal to the 
^-coordinate, one can write ^ as 

V(x, z, t) = V (z, t) + e~^ 4t {^ 1 {z)e- i ^ t - k ^ + <5> 2 (z)e i{uJt ~ k * x) ) , (34) 

where the structure of $i and $2 allows to distinguish between bulk and surface modes. 



The linear variation of the density (transition density), associated with the solution (34), is 
given by 

5p = 2\^ (z)\® + (z)cos(ujt - k x x) (35) 
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with $+ = $1 + $ 2 - 



From the equation of motion (|3T| ) we obtain linear equations for the functions &i(z) 
and $2(2), which, without any loss of generality, can be chosen real. These functions are 
expanded on a basis of eigenstates of the static one-body Hamiltonian Hq in order to get a 
matrix equation which is solved numerically by direct diagonalization. 

In order to write the RPA equations one needs the explicit form of the functional S. The 
functionals used so far in liquid helium have the form |T^Jl^ -[l6| 

So = ^|V*| 2 + VM , (36) 

where V[p] is a velocity independent functional. An accurate phenomenological functional of 
this form [|17| (hereafter named Orsay-Paris functional) has been used, for instance, in Ref. 



[15 1 to study the dispersion of collective modes in helium films. The same functional has 
proven to be quite reliable in the study of the equilibrium configuration of several helium 
systems, such as the free surface [[TT|] , films flIB| , droplets [19|], and vortices |2( 



With functionals of the form ( |36"D the dispersion law predicted by Eq. ( pl|) has the form 

m\x{k)\ 

where \ is the static response function. The quantity x{k) is a key ingredient for calculations 
with density functional theory at zero temperature. A major advantage of liquid ^He is that 
x(k) is well known experimentally |21 . 



It is worth noticing that the dispersion law ([37]) does not reproduce correctly the ex- 
perimental phonon-roton dispersion, especially in the roton region. This is due to the fact 
that functional ( p6|) ignores backflow effects. To understand better this point it is useful to 
rewrite Eq. ( |3TD as the ratio between the energy weighted and the inverse energy weighted 
moments of the dynamic structure function: u 2 = mi (k) /m_i (A;) , where 

f k 2 
mAk) — 1 duj S(k 1 uj) uo = — (38) 
J 2m 

is the f-sum rule, and 
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m_i(Jfe) = J du S{k,u) uj- 1 = --x{k) (39) 
is the compressibility sum rule. The analysis of the spectra of neutron scattering experiments 



21] shows that, due to the u" 1 factor in the integrand, the collective phonon-roton mode 
almost exhausts the compressibility sum rule for all wave lengths up to about 2.2 A -1 . On 
the contrary, it gives only a fraction (~ 1/3) of the energy weighted sum rule, the remaining 
part being exhausted by high energy multiphonon excitations. This explains why functionals 
of the form fl3"ED, yielding the dispersion law fl3"T|), significantly overestimate the roton energy. 
To overcome this difficulty we add a current-current interaction term to the functional 

£ = £ -jJdT> Vj(\v - r'|) p(r)p(r')(v(r) - v(r')) 2 . (40) 

The new term is Galilean invariant and depends on local changes in the fluid velocity field. 
The new dispersion of bulk excitations, obtained by solving Eq. ( j31|) , takes the form 

^ = ^77Tr[l-p(Vj(k)-Vj(0))} , (41) 
m\x(k)\ 

where Vj(k) is the Fourier transform of the current interaction Vj(r) of Eq. (fTO|). Notice 
that the new term in Eq. ( ]4"0|) does not modify the static response function x(k), which is 
completely fixed by Eq. Notice also that the dispersion law fl4lf) is not affected by the new 
term in the k — > hydrodynamic regime, where it coincides with the traditional phonon law 



uj = kj \Jm\x{0)\. This is an important feature ensured by Galilean invariance. 

Equation (f4l|) explicitly shows that the inclusion of the current term in the density func- 
tional can provide a reduction of the phonon-roton contribution to the moment mi, thereby 
improving the description of the dispersion law. The fact that the new functional no longer 
satisfies the f-sum rule (^) points out in a clear way that the time dependent density func- 
tional theory, or equivalently the RPA, does not account for multiphonon excitations which 
provide the remaining fraction of the f-sum rule. A similar separation between collective 
and multiphonon excitations, in the context of linear response theory, has been developed 
by Pines and collaborators (see for example Ref. |2! 
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Let us briefly discuss the criteria used to choose the new functional, whose explicit form 
and parametrization is discussed in Ref . || . The static part V has a form similar to the one 
of the Orsay- Paris functional ]I7|]. It contains a two-body Lennard- Jones interaction and 
a phenomenological density dependent term accounting for short range correlations. The 
density dependent term of the Orsay-Paris functional has been changed in order to reproduce 
more accurately the experimental value of the static response function x(k) in the roton 
region. The current-current interaction Vj is then chosen to reproduce phenomenologically 
the phonon- roton dispersion in bulk liquid |23| through Eq. (fill). Its form is shown in Fig. |2|. 



After fixing the density functional to reproduce the bulk spectrum, one can solve the 
equations of motion in systems with non uniform density. We will focus here on the problem 
of the free surface. We choose, for numerical convenience, a slab geometry (liquid between 



two parallel surfaces) and look for solutions of the form (|34|) . All the states are confined in a 
box larger than the slab thickness, and, consequently, even the continuum of states outside 
the slab is discretized. We take typically a basis of more than 50 eigenstates of H , to cover 
all the relevant part of the spectrum. 

The main features of the spectrum, extrapolated to the semi-infinite system, are shown in 
Fig. [5f The dispersion of bulk phonon-roton modes is the upper solid line; the agreement with 
the experimental data follows from the choice made for Vj. The surface mode corresponds 
to the lowest solid line. An important feature emerging from the present calculation is 
the deviation of the ripplon dispersion from the hydrodynamic law. This is in agreement 
with the experimental data for surface excitations in helium films [|Ij]. The horizontal line 
corresponds to the threshold for rotons having energy close to A = 8.65 K and propagating 
at different angles. The dispersion of the surface mode reaches the roton threshold at about 
1.15 A -1 . The threshold u = |/i 4 | + k 2 x j2m for creating excited states outside the liquid, 
at zero temperature, is also shown. The value at k x = coincides with the experimental 
chemical potential, \/jl4\ = 7.15 K. 

To appreciate the effect of the backflow term included in the functional, the spectrum 
obtained without backflow (Vj = 0) is plotted in Fig. [|. The roton minimum is pushed up 
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significantly, from 8.65 K to about 14 K, as follows from Eq. fl4T|), while the dispersion of 
the surface mode is less affected by backflow. As a consequence the branch of the surface 
mode reaches A near the roton minimum. The spectrum in Fig. is similar to the one of 



previous calculations on slabs p4] , |i~5|| and films [ p~5] , p~3|] . We note however that, even without 
backflow, density functional theory, once the value of x(k) is properly accounted for, gives 
better predictions than the original Feynman approximation. The latter, in fact, gives the 
roton minimum at about 20 K. The reason is that the phonon-roton dispersion in Feynman 
approximation is fixed by the ratio uo = mi(k)/mo(k) between the energy weighted and non 
energy weighted moments of the dynamic structure function. This ratio should be compared 
with prediction ([37]) of density functional theory without backflow. The moment m is 
more affected than the inverse energy weighted moment m_x = — (l/2)x{k) by multiphonon 
excitations. It then follows that the density functional approach gives a significantly lower 
value for the energy of the collective mode (14 K at the roton minimum). 

Due to the finite thickness of the slab, with parallel surface at z = ±L, one can distinguish 
between odd and even states. In Fig. we show a typical spectrum with the excited states 
of a slab of thickness 2L = 50 A. For simplicity, only even states are shown in the figure. 
Each state, at given energy and parallel wave vector, is a stationary superposition of bulk 
excitations, surface excitations and free atoms, with the proper matching conditions. For 
instance, the almost orizontal lines are mainly bulk rotons reflecting at the slab surfaces; the 
gap between these states is fixed by the slab thickness. The oscillations of the roton lines, 
below the maxon region, are peculiar of the slab geometry; they originate from a periodic 
splitting of even and odd states. The period of such oscillations is predicted by the theory 
of roton scattering developed in Sect. to be Ak x ~ pin / (2Lk x ) . The energy of the surface 
mode is instead practically independent of L. Above the ripplon mode there is also an 
excited surface mode, whose dispersion reaches the roton threshold at about k x = 0.7 A . 

The solution of the RPA equations also allows to compute the dynamic structure function, 
defined by 
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S(k,cu) = J2\(n\pl\0)\ 2 5(uj-uj n0 ) (42) 

n 

where p\ = X^-exp^k • r^] is the usual density operator, and k = (k x ,0,k z ). An example 
is given in Fig. |6] where we show the results for k x = 0.8 A -1 and k z = in a slab 50 A 
thick. The ripplon corresponds to the peak at lowest energy. The other peaks are excited 
surface modes, as the one at about 9 K, and bulk modes, the most important one being the 
high energy phonon at about 13 K. The figures reveals that, even for a rather thick slab, the 
contribution of the surface mode to the total strength is significant. It is also interesting to 
explore the dependence of the relative strength on the scattering angle, i.e. as a function of 
k z for a fixed k x . In Fig. [7] the relative contribution of the surface mode to the total strength 
is shown for k x = 0.8 A -1 . One notes that the most favorable condition for exciting ripplons 
is to induce momentum transfer parallel to the surface (k z = 0). 

In Fig. |H the quantity $ + , characterizing the transition density (see Eq. fl35|)), is plotted 
in three cases: a surface mode (a), free atoms coupled to high energy phonons (b), free atoms 
coupled to rotons (c). In Fig. |^ we show in detail what happens to the surface mode when 
its energy approaches the energy of the lowest roton mode. The surface mode transforms 
continously in a bulk mode, and vice-versa. The gap between the two levels tends to vanish 
when the slab thickness increases. The resulting hybridization mechanism then agrees with 
the predictions of the scenario (1) discussed in Section |I| (see Eqs. (|l6j) and (|25|) ), with 
k t ~ 1.15 A" 1 . 

In Sec. m we have shown that rotons close to A, according to Eqs. (^0]) and (|2l|), are 
totally reflected with a change-mode reflection, as in Fig. |]. The roton wave function for 



energy close to A is explicitly written in Eq. (0). In the slab geometry one has to account 



for R + and R~ rotons propagating in both directions and reflecting on two surfaces, at 
z = ±L. The corresponding roton wave function is particularly simple at k x = 0, where it 
takes the form 

7T ' Z 

$! r (z) OC COs(p Q z) COS (7^) • ( 43 ) 

This picture is well confirmed by the numerical solution of the equations of motion. The 
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comparison between the result of the numerical calculation and formula fl43|) is shown in 
Fig. [H]. The two curves are almost indistinguishable, showing that the density functional 
approach properly accounts for the change-mode reflection of rotons on the free surface. 

We note however that result (pip , yielding change- mode reflecton, is not valid when 
a(u),k x ) = 0. We can study the behaviour of rotons also in this case using the arguments of 
Sect. U In fact, for uj > A, the line where a = is fixed by the equation 

-b(k x ) . (44) 



e Q (k x ) - LO 

It starts at the point k t and corresponds to the continuation of the ripplon dispersion above 
A (this continuation is clearly visible in our numerical results of Fig. This line is singular 



for the reflection of rotons at the surface. It is clear from Eqs. ( |19|) and fl20[ ) that, when 
a = 0, one has 

R 1 ~ , R 2 ~ 1 . (45) 

Thus the normal reflection dominates near the singular line. 

We note again that in a slab of finite thickness 2L the roton levels are discretized. For 
the lowest levels near A and far enough from the singular line a = 0, the change-mode 
reflection is dominant. The condition of energy quantization has the form 

2k z L + 28 = wr , n = 1,2,3,... (46) 

where 8 is a phase shift at the reflection, defined by 

R 1 {u) = -\R 1 \exp[2iS\ . (47) 

In the region under consideration R\(u) ~ —1 and one can choose 5 = 0. Approaching the 
singular line a = the quantization condition (f46"D is no longer valid. However, one can 
easily discuss the behaviour of the phase during the crossing of this line. In fact, comparing 
Eq. (|4~7D with Eq. (|19D , one can write 

tan(25) = -— . (48) 
a 
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Crossing the singular line from left to right a changes sign, the change-mode reflection 
amplitude is again equal to —1, but the phase 5 has increased by ir/2. This means that the 
orizontal line &( n )(k x ) for a roton on the left is matched to the line oJ^ n+ i){k x ) on the right of 
the singular line a = 0. These jumps are clearly seen in Fig. [5], just above the point where 
the ripplon branch reaches A. 

We finally note that excited states in the region between the minimum roton energy 
and the maxon energy have a finite amplitude outside the slab, resulting from the coupling 
between rotons and evaporated atoms (see Fig. |8|). A detailed analysis of this coupling can 
be done following the stabilization method recently proposed in Ref. |[25|| . One solves the 
eigenvalue problem in a large box of variable size, containing the slab. The energy of the 
excited states varies with the box size. From the analysis of such energy dependence one can 
extract information about the width of states localized in the slab, i.e., of rotons. A typical 



picture is given in Fig. [II]. The change of slope of each curve corresponds to the coupling 
between free atoms and rotons, related to the phenomenon of quantum evaporation. A 
detailed study of this problem is in progress. 



IV. CONCLUSIONS 

In this work we have provided a detailed discussion of the ripplon dispersion of superfluid 
^ie in the regime of short wave lengths. The main results can be summarized as follows: 

a) We have shown that the hybridization mechanism between ripplons and rotons is 
dominated by the reflection properties of rotons from the surface. We have discussed different 
scenarios for the behaviour of the ripplon dispersion near the threshold uo = A, where A is 
the usual roton gap. 

b) We have calculated numerically the dispersion of ripplons by solving the equations of 
time dependent density functional theory in thick slabs. A phenomenological backflow term 
has been included in the functional in order to provide a realistic dispersion of bulk rotons. 
The resulting ripplon dispersion turns out to be in reasonable agreement with experiments. 
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Our calculations confirm in an explicit way the general reflection mechanism exhibited by 
rotons. 

An extension of the present analysis to investigate higher energy processes associated 
with the phenomenon of quantum evaporation will be the object of a future paper. 
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APPENDIX: 

In this appendix we provide a simple argument to derive expression ( ^3|) for a(k x , u). To 
this purpose one must solve the Schrodinger equation for the roton wave function in the 
domain z <C /t _1 including the boundary. For simplicity, since the wave function \l/ r can be 
chosen real, we take / = fi = f£, so that 

/'(*) 



a(u, k x ) 



m 

The Schrodinger equation for / has the form 



(Al) 

z=0 



(H f -A)f = , (A2) 

where the Hamiltonian Hf(k x ) contains the interaction with the surface. Notice that in 
Eq. ( |A2|) we have put u = A since the cu-dependence of the solution is smooth and can be 
ignored in first approximation. Let the value of a predicted by Eq. ( |A2|) be b(k x ). 



The ripplon dispersion is defined by the Schrodinger equation for the ripplon wave func- 
tion <p 

{H* -w)<f> = (A3) 

whose solutions define the unperturbed ripplon dispersion €o(k x ). 

Hybridization means that on the right hand side of Eqs. ( |A2|) and 
terms linear in and /, respectively: 

(H f - A)f = gC<f> (A4) 



A3) we should add 



(H (t> -uj) ( j> = g£ 1 f (A5) 

where L and C\ are some operators, while g is a coupling constant. Clearly the last equation 
has a pole when uj is equal to the eigenvalue eo(k x ) of H^. Near the pole one can write the 
solution as 

U Pif (A6) 



e (k x ) - uj 

where the quantity £2/ has no dependence on uj. Equation (|A4j) then becomes 

{Hf — A)/ = —npr CC 2 f . (A7) 

e {k x ) - uj 



The solution of Eq. (|A7|) , needed to calculate at(u, k x ), is difficult in general. However it is 
evident a priori that the solution depends on ui only through the combination g 2 / (eo(k x ) —uj). 
One can finally calculate a(u, k x ) by treating the contribution of the right hand side in Eq. 
( |A7| ) as a perturbation. We then obtain the desired result 

Q 2 

a(u, k x ) = b(k x ) + —f- , (A8) 

e {k x ) - uj 

used in Sect. [H|. 

A non perturbative approach to the same result can be also formulated using Green's 
function techniques. 
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FIGURES 

FIG. 1. Schematic picture of the change- mode reflection. A roton with positive group velocity 
(R + ) reaches the free surface and transforms in a roton with negative group velocity (R~)- 

FIG. 2. Effective current-current interaction (adimensional) entering functional (40) 

FIG. 3. Dispersion relation of bulk and free surface excitations. The present result for the 
surface mode (lowest solid line) is compared with the experimental data on films [1] (squares) and 
with the hydrodynamic dispersion of ripplons (dots). The bulk phonon-roton branch (upper solid 
line) is compared with the experimental one (circles) . The threshold for roton states with different 
k x is shown as the horizontal line. The threshold for the emission of atoms into the vacuum is also 
shown (dot-dashed curve). 

FIG. 4. Same as in Fig. 2 but with the bulk and surface mode dispersions calculated without 
the current-current interaction in the density functional (Vj = 0). 

FIG. 5. Dispersion relation of the excited states of a slab 50 A thick, as a function of the 
parallel wave vector. 

FIG. 6. Dynamic structure function (in arbitrary units) for k z = and k x = 0.8 A -1 in a slab 
50 A thick. The ripplon mode corresponds to the peak at lowest energy. 

FIG. 7. Ratio between the strength of the surface mode and the total strength at k x = 0.8 A -1 
as a function of k z . 

FIG. 8. Three examples of excited states of a slab 50 A thick. The dashed line is the static 
density profile. The solid line corresponds to the numerical solution of the equations of motion, 
at given values of energy and parallel wavevector. The three states have the same k x but different 
energy: a) the lowest energy surface mode; b) free atoms outside the slab coupled to phonons; c) 
free atoms outside the slab coupled to rotons. 
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FIG. 9. Detail of the repulsion between the surface mode and the bulk roton mode near the 
threshold energy A, for a slab 50 A thick. The form of the lowest eigenstate at two different values 
of k x is also shown in the upper part of the figure. 

FIG. 10. Wave function of the first excited state above the roton threshold at k x = 0. It 
corresponds to a superposition of R + and R~ rotons which undergo a change-mode reflection at 
the surface. Solid line: numerical solution of the equations of motion. Dashed line: analytic 
expression (43). Short dashed line: static density profile. 

FIG. 11. Energy of the excited states at k x = 0.4 A -1 as a function of the size of the compu- 
tational box, containing a slab of fixed thickness (50 A). The almost horizontal part of each curve 
corresponds to roton-like states inside the slab; the part with negative slope (proportional to L7 2 ^) 
corresponds to free atom states outside the slab. The form of the level repulsion is connected to 
the probability for a roton to decay by ejecting a free atom from the liquid. The energy of the 
surface modes (not shown in the figure) is practically independent on L^ Qyi . 
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